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Abstract 

The search for lepton-flavor-violation process e + e~ — ► e/i in the energy region ^/s = 984 - 1060 
MeV with SND detector at VEPP-2M e+e" collider is reported. The model independent 90% CL 
upper limits on the e + e~ — > e/i cross section, <r eA1 < 11 pb, as well as on the corresponding <j> — > e/x 
branching fraction, B(cf> — ► e/i) < 2 x 10 -6 , for the final particles polar angles 55° < 9 < 125°, were 
obtained. 
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For the most of fundamental fermions (quarks and neutrinos) the processes with flavor 
violation, quarks decays and neutrinos oscillation, are known. At the same time the LFV 
processes with charged leptons has never been observed. Theoretically the processes of this 
kind are not strictly forbidden and can occur in many extensions of the Standard Model. 

For the LFV hunting, the decays of \i and r leptons, as well as of the Z-boson and 
of various quark-antiquark mesons (K, B, D,rj, J/ip,T), along with a conversion process 
fj,N — > eN are used [l[ Q]. The annihilation processes e + e~ — ► e/i 7 er, /it are also suitable 
for this purpose. Theoretically these processes and related gauge boson and vector meson 



decays were studied, for example, in [3J. On the experimental side, the searches for the 
decays J/ip — > e/i, er, fir Q] , T — > fir j^, Z — > e/i, er, /ir jf|, as well as for the annihilation 
processes e + e~ — ► er, /ir in the T(4S r ) energy domain J7|, and for the processes e + e~ — > e/i, 
er, fir in the energy region ^Js = 189 - 209 GeV [8j were performed. However, in the 
energy region below the J/ip production threshold such studies were not done yet. In the 
0(lO2O)-meson energy domain, it is possible to search for the LFV process e + e~ — > e/i and 
the corresponding decay (p — > e// (FigJTJ) . 




FIG. 1: The diagrams of the e + e — ► e[i process. 

Existing stringent bounds on LFV /x — > 3e decay can be transformed to a severe constraint 
on the two-body <fi —>■ e\x branching fraction: B(<p — > e/i) < 4 x 10 -17 unless some magic 
cancellations take place in the \i — > 3e decay amplitude [9j]. At first sight, such a strong 
constraint makes doubtful any experimental effort to search this decay. However, the magic 
cancellations mentioned above, although unlikely, cannot be absolutely excluded. 

This work reports the results of studies of the process e + e~ — > e/i in the energy region 
v/i ~ 1 GeV with SND detector at e+e" collider VEPP-2M. 
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The SND detector [10] operated from 1995 to 2000 at the VEPP-2M [U| collider in 
the energy range \fs from 360 to 1400 MeV. The detector contains several subsystems. 
The tracking system includes two cylindrical drift chambers. The three-layer spherical 
electromagnetic calorimeter is based on Nal(Tl) crystals. The muon system consists of 
plastic scintillation counters and two layers of streamer tubes. The calorimeter energy 
and angular resolutions depend on the photon energy as <te/E(%) = 4.2%/ y/E(GeV) and 
°<j>,e = 0.82° /y^E(GeV) © 0.63°. The tracking system angular resolution is about 0.5° and 
2° for azimuthal and polar angles respectively. 

This work is based on the data collected in the scans of the 0-meson energy region. The 
total integrated luminosity used is IL = 8.5 pb _1 . The luminosity was measured using the 
process e + e~ — ► e + e~ with the accuracy of about 2%. 

In the reaction e + e~ — > e/i the final particles are detected by the tracking system and 
have substantively different energy depositions in the calorimeter. The muon system detects 
muons with a probability of greater than 90%, while electrons are detected by this system 
with the probability of less than 0.2%. To search for e + e~ — > e/i process, the so called 
collinear events containing two charged particles were used. We assume that the charged 
particle with higher energy deposition in the calorimeter (particle number one) is an electron, 
while the particle with lower energy deposition (particle number two) is a muon. The events 
were selected using the following criteria (subscripts 1 and 2 denote the particle number): 

1- N c ha = 2, where N c h a is the number of the charged particles originated from the 
interaction point: \zi^\ < 10 cm and ri,2 < 1 cm. Here z is the coordinate of the 
charged particle production point along the beam axis (the longitudinal size of the 
interaction region a z about 2.5 cm), r is the distance between the charged particle 
track and the beam axis in the r — <j) plane; 

2. |A0| = 1 180° - (0i + e 2 )\ < 20°, where 9 is the particle polar angle; 

3. |A0| = 1 180° — |0i — 02 1 1 < 5°, where is the particle azimuthal angle; 

4. 55° < 1>2 < 125°; 

5. the angular region 240° < 0^2 < 300° not covered with the muon system was excluded; 

6. the muon system was hited by the second particle and was not hited by the first one; 
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7. 20 < E{ < 50 MeV, 40 < E^ 1 < 80 MeV and 50 < E{ n < 90 MeV, where E\ are 
the energy depositions in the calorimeter layers, i denotes the particle number and 
j = I, II, III is the layer number; 

8. E[ > 70 MeV, E[ T > 130 MeV and 20 < E[ u < 100 MeV. 

As a result 146 events were selected. The visible cross section (the events number divided 
by the integrated luminosity) varies weakly with beam energy. No contribution from the 4>- 
meson decays — > K + K~ , KgK^, 7r + 7r~7r° is seen. This agrees with the expectations from 
the Monte-Carlo (MC) simulation. The events from the background process e + e~ — > 7r + 7r~ 
can pass the selection if one of the pions looses its energy due to the ionization, while the 
other pion - due to nuclear interactions. 




FIG. 2: The first particle E*/E e for selected events (dots with errors). Solid curve - fit by a sum 
of the distribution for electrons (histogram) and Gaussian. Dashed curve - fit by a sum of the 
distribution for electrons (histogram) and third-order polynomial. 

In order to obtain the cross section of the process e + e~ — > efi in the whole energy region 
y/s = 984 - 1060 MeV, the E*/E e distribution (FigJ2J) was analyzed. Here E* and E e are 
the electron energies measured by the calorimeter and expected from the process kinematics 
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respectively. To obtain the number of e + e~ — > e/i events (N efi ), the E*/E e spectrum was fit 
by a sum of distributions for electrons and background. The distribution for electrons was 
obtained using experimental data. The background was approximated either by Gaussian 
function or by a third-order polynomial. The coefficients of the background function and 
N efl were free parameters of the fit. When the background was approximated with Gaussian, 
it was found that 



N efl = 12 ±\* . 



This corresponds to the upper limit 



iV eM <30 CL=90%. 



In the case of the third-order polynomial it was obtained 



N efM = 7±l 1 . 



The corresponding upper limit is 



N efl < 21 CL=90%. 



The higher limit N efl < 30 was used for the further considerations. 

Tracking system detection efficiency for the e + e~ — > e/i events, E tra ck (cuts 1-5), was 



obtained from MC simulation 
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121 ]. The MC events were generated with 1 + cos 9 
distribution. The detection efficiency obtained for the angular region 55° < 9 < 125° 
actually does not depend on the model of the ^-distribution. It's equal to Etrack = 0.59. The 
experimental and simulated 9, A9 and A<p distributions for the processes e + e~ — * e + e~, 
e + e~ — > 7r + 7r~, e + e~ — > /i + /t~ are in a good agreement [3, Q. The systematic uncertainty 
associated with e track determination is estimated to be less than 3 %. 

The efficiencies of muon and electron detection by the muon system were obtained using 
e + e" — > ji + n~ and e + e~ — > e + e~ data events. The e + e~ — > events were selected 

according to the criteria 1-5 described above. The additional cut 7*1,2 < 0.1 cm was used 
for suppression of the cosmic ray background. The cut 7 was imposed on both particles. 
One particle was required to hit the muon system, while the other particle was used to 
determine the detection efficiency. The residual cosmic background was subtracted using the 
distribution of the ^-coordinate of the particles production point [13J]. Detection efficiency 
£ muon depends on the muon energy. Its value varies from 0.90 to 0.95 with the average value 



of £^ uon being equal to 0.94. The e + e — > e + e events were selected by the cuts 1- 5 and 
the condition E lj2 /E > 0.7. It was found that 1 - e e muon = 0.998. 

Probabilities for muons and electrons to pass condition on the energy deposition in the 
calorimeter were obtained in a similar way. The e + e~ — > fi + fi~ events were selected using 
cuts 1-5 and additional requirements 7*1,2 < 0.2 cm, Ei t2 /E < 0.6. It was required that 
the muon system was hit by the both particles. The cosmic background was suppressed by 
the restriction \t± — r 2 | < 5 ns, where Ti )2 are time intervals between the signals from the 
scintillation counters and the beam collision moment. The probability for muons to pass 
cut 7 was found to be e^ al = 0.86. The e + e~ — > e + e~ events were selected using criteria 1 
- 5 and requiring that the muon system was not fired by any particle and that the energy 
deposition of a randomly chosen particle was greater than 0.85 x Eq. Other particle was 
used to obtain probability to satisfy the criterion 8, el al = 0.70. The values of and e e cal 
do not depend on the beam energy. 




FIG. 3: The visible cross section obtained after additional cut 0.9 < E*/E e < 1.1. The solid curve 
is expected resonance line shape corresponding to the upper limit on B(4> ^ e/j,), dashed curve is 
background approximation. 
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The detection efficiency of the e + e — > e/i process was calculated as follows 



£ efi — Ztrack x £ rnuon X (1 £ muon) X £ cal X 6 



cat 



and the e + e — > e/i cross section as 



cr P . 



ILe efJi 



where N efl < 30, e e ^ = 0.31 (the average value s^uon = 0.94 was used), IL = 8.5 pb l . The 
following upper limit for the angular region 55° < 9 < 125° was obtained 

cr eM < 11 pb CL=90%. 

The upper limit on the <fi — > efi decay was obtained assuming absence of any non-resonance 
contribution and by using the additional cut 0.9 < E*/E < 1.1 (then = 0.64). The 
energy dependence of the visible cross section is shown in Figj3j It was fit by the function: 

2 



(J = £ X (1 + 5 rad ) X 



3,s 



■V '5(0 -»• e+e-)B(<p -f e/i)- 



where (1 + <5 ra d) is the radiative correction factor [jjj], ^(s) is a second-order polynomial 
describing the background, m^, are the 0-meson mass and total width, respectively. The 
branching ratio and the coefficients of the ^(s) were free fit parameters. For the angular 
region 55° < 9 < 125°, it was obtained 

B(<p -> e/i) = (0.0 ± 1.5) x 10~ 6 , 

which corresponds to the upper limit 

B((f) -> e/i) < 2 x 10~ 6 CL=90%. 

The presented upper limits do not depend on the angular distribution of the process e + e~ — ► 
e/i. 
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